Identifying the odd-frequency superconducting state by a field-induced Josephson effect 
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Superconducting order parameters that are odd under exchange of time-coordinates of the electrons consti- 
tuting a Cooper-pair, are potentially of great importance both conceptually and technologically. Recent exper- 
iments report that such an odd-frequency superconducting bulk state may be realized in certain heavy-fermion 
compounds. While the Josephson current normally only flows between superconductors with the same symme- 
tries with respect to frequency, we demonstrate that an exchange field may induce a current between diffusive 
even- and odd-frequency superconductors. This suggests a way to identify the possible existence of bulk odd- 
frequency superconductors. 

PACS numbers: 74.20.Rp,74.25.Fy,74.45.+c,74.50.+r 
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I. INTRODUCTION 

The prevalent symmetry in known superconductors may 
be described as odd under exchange of spin coordinates, and 
even under an exchange of spatial coordinates or an exchange 
of time coordinates of the electrons constituting the Cooper 
pair. The latter condition corresponds to a sign change of 
frequency, after Fourier-transforming the time variables to 
a frequency representation. This symmetry may compactly 
be expressed as an even-frequency singlet even-parity super- 
conducting state (hereafter referred to as the even-frequency 
state). 

However, other types of pairing are also permitted. Among 
these is the so-called odd-frequency pairing state 1 . Such a 
state is potentially of great importance, both from a conceptual 
as well as a technological point of view. From a conceptual 
point of view phase transitions involving pairing of fermions 
form centerpieces of physics in such widely disparate sub- 
disciplines as cosmology, astrophysics, physics of condensed 
matter, physics of extremely dilute ultra-cold atomic gases, 
and physics of extremely compressed quantum liquids. Ex- 
tending the possible pairing states compatible with the Pauli- 
principle will likely have impact on all these disciplines, and 
hence deserves attention. From a technological point of view, 
an odd-frequency pairing state would be a candidate for a ro- 
bust superconducting pairing state capable of coexisting with 
ferromagnetism 2 . Such a material would be extremely im- 
portant, since it would combine two functionalities of major 
importance, namely magnetism and superconductivity. 

The hope of experimentally detecting this type of pair- 
ing was raised by the prediction that an odd-frequency 
triplet even-parity (hereafter simply denoted odd-frequency) 
pairing should be induced by the proximity effect in dif- 
fusive ferromagnet/conventional superconductor junctions 3 . 
Recently, odd-frequency states due to the proximity effect 
have been confirmed in experiments of diffusive ferromag- 
net/conventional superconductor junctions 4 . To date, no bulk 
odd-frequency superconductor has been unambigously identi- 
fied. 

Currently, the heavy-fermion compounds CeCu2Si2, and 
CeRhlns seem to be the most promising candidates for the 



5odd/evcn( X 




<7odd/ even (+x) 



FIG. 1 : A diffusive metal which may be either normal or ferromag- 
netic (with exchange field h) of width d separates two superconduc- 
tors with a phase difference 2\, thus constituting a Josephson current 
(J is the current source). The individual superconductors have either 
an even-frequency or odd frequency bulk symmetry. The exchange 
field may have any orientation in the yz-plane, characterized by the 
angle cj>. 



realization of a bulk odd-frequency state 5 . It has also recently 
been argued 6 that this type of pairing could be realized in 
Naj;Co02, motivated by band-structure calculations and the 
robustness of its superconductivity against impurities. How- 
ever, the experimental reports on the Knight-shift data have 
shown evidence of both singlet 7 and triplet 8 pairing. To re- 
solve the pairing issue, it would be desirable to make clear-cut 
theoretical predictions for experimentally measurable quanti- 
ties that may distinguish the odd-frequency symmetry from 
the conventional even-frequency symmetry. Motivated by 
this, the conductance of a diffusive normal metal (N) /odd- 
frequency junction was recently studied 9 . 

It is well-known that the Josephson current between super- 
conductors with different symmetries is in general inhibited 10 . 
However, it was very recently shown that in the clean limit, 
such a Josephson coupling may be established between such 
superconductors by means of surface-induced pairing com- 
ponents of different symmetry than the bulk state in each 
superconductor 11 . In the dirty limit, it has been demonstrated 
that in the case of an s-wave even- or odd-frequency bulk su- 
perconductor, the proximity-induced pairing component in the 
N will have the same symmetry 12 . These results are valid in 
the absence of an exchange field. Due to the above mentioned 
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properties, one would consequently expect that both the even- 
and odd-frequency symmetries are induced in the normal part 
of an even-frequency /N/odd-frequency junction. While this is 
true, the Josephson current is found to vanish in such a setup. 
If the N is replaced with a diffusive ferromagnet (F), the pair- 
ing components in the F will have the same symmetries as in 
the N case. Surprisingly, the Josephson current is not absent 
in this case. 

In this paper, we report that a Josephson current may flow 
in a diffusive even-frequency /F/odd-frequency junction, and 
argue that this should serve as a smoking gun to reveal the odd- 
frequency symmetry in a bulk superconductor. We also study 
the dependence of the Josephson current on the temperature 
and width of the F, and show that 0-ir transitions take place. 
In the following, we will use boldface notation for 3-vectors, 
. 7 . for 8x8 matrices, . ? . for 4x4 matrices, and for 2x2 
matrices. 



II. THEORETICAL FRAMEWORK 

In order to address this problem, we employ the quasiclassi- 
cal theory of superconductivity using the Keldysh formalism. 
The quasiclassical Green's functions may be divided into an 
advanced (A), retarded (R), and Keldysh (K) component, each 
of which has a 4 x 4 matrix structure in the combined particle- 
hole and spin space. In thermal equilbrium, it suffices to con- 
sider the retarded component since the advanced component is 
obtained by g A = —(p3g R p3)' i , while the Keldysh component 
is given by g K — (g R — g A ) tanh(/fe/2), where j3 is inverse 
temperature. Pauli-matrices in particle-hole x spin (Nambu) 
space are denoted as pi, while Pauli-matrices in spin-space are 
written as r,, i = 1, 2, 3. As shown in Ref^, the symmetry 
of the anomalous Green's function induced in the N through 
the proximity effect by an odd-frequency superconductor is 
also odd-frequency, and similarly for an even-frequency su- 
perconductor. We may write the retarded Green's function in 
the superconductors as 
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for an odd-frequency triplet even parity s-wave symmetry / 
even-frequency singlet even parity s-wave symmetry. Above, 
we have defined 

c = cosh(0), s = sinh(0), 9 = atanh[A(e)/e], (2) 

where e is the quasiparticle energy measured from Fermi level 
while x is the superconducting phase associated with the bro- 
ken U(l) symmetry. The ±-sign refers to the right and left 
side of the F (Fig. [T). The retarded Green's function in the F 
satisfies the Usadel equation^ 

D\7(g R \7g R )+i[er 3 + M,g R } = 0, 
( cos </>T3 + sin (f>T2 



M = h 



V 







(COS (j)T3 + Sin (j)T2 ) 
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where D is the diffusion constant, h is the exchange field, and 
T denotes the matrix transpose. An analytical solution of this 



equation is permissable when it may be linearized, which cor- 
responds to a weak proximity effect. This may be obtained in 
two limiting cases: 1) if the barriers have low transparency or 
2) if the transmission is perfect (ideal interfaces) and the tem- 
perature in the superconducting reservoir is close to T c , such 
that | A | is small. We will here mainly focus on the low trans- 
parency case 14 since it might be difficult to experimentally re- 
alize small and highly transparent junctions for observing the 
Josephson current. In this case, we make use of the standard 
Kupryianov-Lukichev boundary conditions^ 



2jd(gd x g) = ±[g,g R{L) {±x)}\ x =d,o, 



(4) 



where 7 is a measure of the barrier strength and <7r(d denotes 
the Green's function in the superconductor on the right (x — 
d) or left (x = 0) side of the F, respectively. Defining the 
vector anomalous Green's function 

f =[/|-/T,-i(/T+/4)> 2 /t]/2. (5) 
and the matrix anomalous Green's function in spin-space, 

z -</•+'■***- Gt* //'//■)' ( » 

the linearized Green's function in the F reads in total 
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Note that Eq. (0 contains both equal-spin and opposite-spin 
pairing triplet components in general. In the special cases of 
= and = 7T, the equal-spin pairing components f a 
(o - =Tii) vanish. 



m. RESULTS 

We now provide the analytical results for the Josephson cur- 
rent in an even-frequency/F/odd-frequency junction. To begin 
with, we will consider the case <f> = to emphasize our main 
result: namely a field-induced Josephson effect. We will then 
investigate the effect of a change in orientation of the mag- 
netization in the ferromagnetic layer. Changing the orienta- 
tion, i.e. (f>, in an even-frequency/F/even-frequency junction 
has no effect since the order parameters in the superconduc- 
tors in that case are isotropic. However, for an odd-frequency 
triplet superconductor the order parameter has a direction in 
spin-space. We here consider an opposite-spin pairing 5*^ = 
order parameter without loss of generality, and then proceed 
to vary the orientation of the exchange field. 



A. Field-induced Josephson current 

For = 0, one readily finds that f a = 0. The linearized 
version of equation Eq. (O is then given by 



Dd 2 J ± + 2i(e±h)f ± =0, 



(8) 
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with the general solution 

f ± (x) = A ± e lk±x + B±e~ lk±x , 



k± = y/2i(e ± h)/D, 



(9) 



with the definition f± = /t±/ s . Employing the boundary con- 
ditions, one finds that the anomalous Green's function reads 

f± = ^ (CL + P± , )T ; SLe " X + B ± e-**', (10) 
(p± - c L )e~ lk ± x 

where we have introduced the auxiliary quantity 

B ± = [ SR c« ± s L e< k ± d -^(p± + cr)/( P± - c L )] 
x [e lk±d (p± + c L )( j o± + c R )/( j o± - cl) 



e- lk±d {p±-c R )X 



(11) 



Moreover, we have introduced p± = r )d\k± while the sub- 
scripts 'L' and 'R' denote the left and right superconductor 
for the coefficients c and s. Once g R has been obtained, the 
Josephson current may be calculated by the formula 

j(a:) = -(N F eDx/4) J deTr{p 3 (gd x g) K } 

/CO 
deRe{M + (e) +M_(e)} 
-OC 

x tanh(/3e/2), (12) 

with the definition 

M±(e) = [f±(-e)]*d x fr(e) - f±{e)d x [f^{~e)]*. (13) 

The normalized current density is defined as 

Hx)/Io = Mi(x,x)\/(N F eDA 2 ), (14) 

which is independent of x for x £ [Q,d\, and the critical cur- 
rent is given by I c = 

At this point, we are in a position to compare the results 
for the even- and odd-frequency case against each other, to 
investigate how the different symmetry properties alter the 
Josephson current. Similarly to Ref. 12 , we will model the odd- 
frequency gap by 



A(e)= e /[l + ( e /A n 



where n = 2,4,8,.., which exhibits the low-energy be- 
haviour considered in Ref. 9 . We choose specifically n = 2, 
and underline that none of our qualitative conclusions are al- 
tered by choosing n = 4, 8. It is natural to begin look- 
ing for signatures of the odd-frequency symmetry as probed 
by the Josephson current in the simplest case of an even- 
frequency/N/odd-frequency junction. The proximity-induced 
anomalous Green's function in the N will have a contribution 
from both the even- and odd-frequency symmetry!-. How- 
ever, we find that no Josephson current may flow in such a 
setup. 



This may be understood by considering the boundary 
conditions at each interface (see Fig. |2). At the even- 
frequency /N interface, the odd-frequency triplet component of 
the proximity-induced Green's function is absent since pen- 
etration into the even-frequency superconductor is prohib- 
ited. Similarly, at the N/odd-frequency interface, the even- 
frequency singlet component of the Green's function vanishes 
for the same reason. 



Therefore, the Josephson current is zero in this type of junc- 
tions since the current-carrying Green's function in the N in- 
duced by the left superconductor is absent at the right super- 
conductor, and vice versa. This is a direct result of the differ- 
ent symmetries of the even- and odd-frequency superconduc- 
tors. The Josephson current can be divided into the individual 
contributions from / s and f t (cross terms vanishes), and for 
each component the coherence is lost since the even- and odd- 
frequency pairings cannot reach the opposite interface. 



Analytically, one can confirm from Eq. ( fTOb and ( fTTT ) 
that the Green's functions providing the critical current (at 
X = vr/4) satisfy f±(e) = e) for any x in the 

even-frequency/N/odd-frequency case, which upon insertion 
in Eq. ( fl~3l ) yields j(x) — 0. In even-frequency/N/even- 
frequency and odd-frequency/N/odd-frequency junctions, one 
can in a similar manner confirm that A(e) = A*(— e)e~ lkd 
and B(e) = B*{— e)e lk , which leads to a finite value of the 
Josephson current. This may be seen for instance at a; = d/2 
upon substitution into Eq. ( fT3l ), where f±(e) = s) and 



(15) d x f±(e) =-B x m- 



■e . 



If we instead consider an even-frequency/F/odd-frequency 
junction, the proximity-induced anomalous Green's function 
in the F will have the same symmetries as in the N case. Re- 
markably, we find that in this case, however, a Josephson cur- 
rent is allowed to flow through the system. The reason for 
why a Josephson current is present when the N is replaced 
with a F, is that the exchange field allows for both the singlet 



and triplet components to be induced throughout the F due 
to the mixing of singlet and triplet components by breaking 
the symmetry in spin space, regardless of the internal sym- 
metry of the superconductors 1617 as shown in Fig. [2] As 
a result, we suggest that as a way to unambigously identify 
an odd-frequency superconductor, a coupling through a N to 
an ordinary even-frequency superconductor should not yield 
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FIG. 2: (Color online) Plot of the singlet / s (thick, blue lines) and triplet f t (thin, red lines) anomalous Green's function induced in a N (upper 
panels) and F (lower panel) in proximity with an even-frequency and/or odd-frequency superconductors, using \ = 7r /4 and e/Ao = 0.1. 
The solid lines are the real parts of the anomalous Green's functions, while the dashed lines are the imaginary parts. To ensure that / S(t <C 1, 
we take d/£ = 50 in the N case and d/£ = 2, h/Ao = 50 in the F case. Here, £ = ■^/D/(2wT c ) is the superconducting coherence length. 
While a Josephson current flows when the symmetries of the superconductors are equal in the N case, it vanishes when the superconductors 
have different symmetries. This is because the triplet and singlet components are absent at the interface of the even- and odd-frequency 
superconductor, respectively. Replacing the N with a F establishes a Josephson coupling between superconductors with different symmetries. 



any Josephson current while replacing the N with an F should 
allow for the current to flow. One could argue that this is 
precisely the case also for even-frequency triplet odd-parity 
superconductors 17 . This symmetry is nevertheless easily dis- 
tinguished from an odd-frequency s-wave symmetry since the 
former is highly sensitive to impurities while the latter does 
not suffer from this drawback. 

We have investigated the behaviour of the critical current 
when coupling an even- and odd-frequency superconductor 
through a F, and the results are shown in Fig. [3] for a repre- 
sentative choice of parameters. To make the numerical calcu- 
lations stable, we have added a small imaginary number— in 
the quasiparticle energy, e — > e + id with S — O.OlAo, and 
we fix 7 = 5 (also for Fig. |2|. One observes the well-known 
0-7T transitions 19 upon increasing with the width d of the F 
[Fig. [3^)]. Similarly, we also find that Q-tt oscillations occur 
as a function of temperature. The current-phase relationship 
is sinusoidal as usual in the linearized treatment. The width 
of the junction is measured in units of the superconducting 
coherence length 25 £ = y/D/(2irT c ). 

As a further probe of the odd-frequency symmetry in a 
bulk superconductor, we investigate the local density of states 
(LDOS) in the F. The LDOS is altered from its normal state 
value due to the proximity-induced anomalous Green's func- 
tion in the F. Consider Fig. for a plot of the deviation 



SN from the normal state LDOS in an even-frequency/F/odd- 
frequency junction. The deviation is given by the formula 

SN(e) = (MVl + f*(e)[f-*{-e)]*} ~ ( 16 > 

cr = ± 

under the assumption of a weak proximity effect. As is seen, 
the LDOS is enhanced at e = due to the presence of 
odd-frequency pairs 12,20 , and the usual peak arises at e = 
Ao. While the odd-frequency/F/odd-frequency case exhibits 
the first property, and the even-frequency/F/even-frequency 
case the latter, the even-frequency/F/odd-frequency junction 
is characterized by the fact that both of these features appear 
in SN. This could serve as an identifier of the odd-frequency 
symmetry in conjunction with the other properties we have 
analyzed here. 

It is interesting to observe the field-dependence of the crit- 
ical current in an even-frequency/F/odd-frequency junction, 
shown in Fig. [4] To gain access to the regime of a very weak 
or absent exchange field, we must choose the width d suffi- 
ciently large to ensure a weak proximity effect. In Fig. @] 
we plot the current as a function of h for several values of 
d/£. One observes that the critical current goes exactly to 
zero at h = 0, while a current is induced for non-zero val- 
ues of h. A maximum peak appears for very weak exchange 
fields, and the critical current oscillates with increasing field- 
strength. Although the short-junction regime is not accessible 
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FIG. 3 : (Color online) a) The critical current (% = 7r/4) as a function 
of d/£ for T/T c = 0.001. b) Deviation from the LDOS at x = d/2 
for an even-frequency/F/odd-frequency junction with \ = t/4 and 
h/ Ao = 50. The zero-energy peak is clearly discernible, originating 
with the odd-frequency pairs. 



for very weak exchange fields due to the linearized treatment, 
it follows from our analytical expressions that the current is 
absent for any choice of d as long as h = 0. 



B. Case with 0/0 

We now proceed to consider the effect of rotating the ex- 
change field in the ferromagnetic region. In effect, we allow 
for cj) 0. Now, the equal-spin pairing components f a are 
in general non-zero, and we find the following four coupled, 



d/£ = 30 
d/£ = 40 
d/£ = 50 




FIG. 4: (Color online) Field dependence of the critical current for 
several values of d/£ with T/T c = 0.001. 



linearized Usadel equations: 

Dd 2 J ± + 2i(e + ftcos0)/± ± Jisin0(/ T + /J = 0, 

Dflg/, + 2ie/ CT - 2hshx4>f s = 0, a =t,| • (17) 

The general solution for the anomalous Green's functions may 
be written compactly as: 

/±(a;) = A(x)(ai + a ; ) + B(a:)(6 T + 6 X ) 

+ c 5 R±{x) + c 6 P±(x) + c 7 R±(x) + c 8 P±(x), 

f a = a a cos Lqx + b a sin LqX + C5 (1 sin L + ir sin </>) 
+ C6(icosZ/+a;sin0) + cj{~ isinL-xsmcj)) 
+ cs(— lcosL-ssin^), <7=t,J,. (18) 

In the above, {a CT , C5, eg, C7, eg} are unknown coefficients 
to be determined from the boundary conditions in the prob- 
lem. Also, we have introduced the following auxiliary quanti- 
ties: 



A(x) 
R%(x) 

L 



1 cos Lqx sin ( 



B(x) 



isinio a; sin( 



2 cos 4> 2 cos < 

sinL-|-a;(±cos0 + a), a = ±1 
cos L±x(± cos (p + a), a = ±1 



(1 + , l± = (i±i; 



(19) 



When the left superconductor has an even-frequency symme- 
try and the right superconductor has an odd-frequency sym- 
metry, the boundary conditions yield at x = 0: 



2jdd x f ± = 2/±c L T2 SL e- lx , 
2ldd x f„ = 2f a c Ll 



(20) 
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while at x = d one obtains 

2 7 dd x f± = 2s R e^ - 2f±c R , 

2^dd x f a = -2/ ct cr. (21) 

The presence of equal-spin pairing components slightly mod- 
ified the expressions Eq. (fT2l and Eq. ( fT6b for the density of 
states and the Josephson current, respectively. We now obtain 

8N(s) = Y^Wa(s) l]/2, 

cr 

N a (s) = Re{{l + U{e)f* a {-e) + [f t (e) 

+ af s (e)][f:(~e)-af:(-e)}}^ 2 } (22) 

for the density of states, while the Josephson current is calcu- 
lated according to 

/oo 
deRe{M+(e) + M_(e) 
-OO 

M T (e) + M;(e)} x tanh(/3e/2), (23) 
with the definition (cr =T> 4) 

M a (e) = [fA-e)}*d x Me) - f a (e)d x [f a (~e)}* . (24) 

Let us first address the issue of how the zero-energy peak 
in the DOS treated earlier is affected by a rotation of the ex- 
change field. In Fig. [5] we plot the deviation SN(0) from the 
normal-state zero-energy DOS as a function of the misorien- 
tation (j) for d/£ = 2 and h/Ao = 30. As </> is increased, it is 
seen that SN(0) grows rapidly. As it becomes comparable to 
the normal-state DOS in magnitude, the linearized treatment 
of the Usadel equations becomes less accurate, denoted by the 
symbols in Fig. [5] Nevertheless, the trend seems clear: the 
zero-energy DOS reaches a maximum at <fi = n/2. Also, we 
have demonstrated (not shown) that the characteristics in Fig. 
|5J)) remain the same for all <f>. In particular, the zero-energy 
peak is not destroyed by changing cf>. 

Examining the magnitude of the anomalous Green's func- 
tion numerically, we find that the weak-proximity effect as- 
sumption \f\ < 1 becomes poor for energies close to zero 
when <j> is close to tt/2. Therefore, this parameter regime is 
strictly speaking inaccessible within our linearized treatment. 
Note that no such problem occurs when <fi = 0. However, we 
will assume that the linearized treatment is still qualitatively 
correct when <j> is close to n/2 in order to investigate how 
the critical current depends on the orientation of the exchange 
field orientation. 

In Fig. [6] we plot the variation of the critical Josephson 
current as a function of the orientation c6 of the exchange 
field. If the two superconductors have a conventional even- 
frequency symmetry, the Josephson current is completely in- 
sensitive to the orientation of the exchange field. This is rea- 
sonable, since the superconducting order parameter in this 
case is spin-singlet and has no orientation in spin space. Note 
that a magnetic flux threading a Josephson junction in gen- 
eral gives rise to a Fraunhofer modulation of the current as 
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FIG. 5: (Color online) Plot of the deviation from the normal-state 
zero-energy density of states in the middle of the ferromagnetic re- 
gion (x/d — 0.5). The symbols denote approximately the region 
where the linearized treatment becomes less accurate, i.e. the anoma- 
lous Green's functions / no longer satisfy |/| -C 1. 

a function of the flux. We here neglect this modification by 
assuming that the flux constituted by the ferromagnetic region 
is sufficiently weak compared to the elementary flux quantum. 
This is the case for either a small enough surface area or weak 
enough magnetization (the energy exchange splitting may still 
be significant). 




0/7T 

FIG. 6: (Color online) Plot of the modulation of the Josephson cur- 
rent with orientation <j> of the exchange field in the ferromagnet. The 
symbols correspond to an even-frequency/F/odd-frequency junction, 
while the dashed line denote an even-frequency/F/even-frequency 
junction with h/Ao = 30 and d/£ = 2.0. 

The situation is quite different when one of the supercon- 
ductors has an odd-frequency symmetry. In this case, the 
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Josephson current is sensitive to the orientation 4> of the ex- 
change field, and displays the behaviour shown in Fig. [6] 
The reason for this is that the order parameter in the odd- 
frequency triplet case has a direction in spin space, here cho- 
sen as opposite-spin pairing (along the z-axis). For <j> = 0, 
the Cooper pairs are opposite spin-paired relative to the ex- 
change field, while for = tt/2 the Cooper pairs are equal- 
spin paired relative to the exchange field. Tuning the rela- 
tive orientation of the exchange field and the superconduct- 
ing order parameter in the odd-frequency superconductor is 
thus seen to lead to the possibility of controlling the magni- 
tude of the Josephson current. In an experimental situation, 
only the orientation of the exchange field is probably alter- 
able. In Fig. [7] we show the current-phase relationship for an 
even-frequency/F/odd-frequency Josephson junction for sev- 
eral values of 4> to show that although it remains sinusoidal, 
it is shifted from ~ sin0 to ~ sin[0 + a (</))] where a(4>) is 
nonzero for h ^ 0. 

Interestingly, for cf> = it/ 2 the current vanishes completely 
as seen in Fig. [6] but it is nonzero otherwise. This can be 
understood by studying the contribution to the Josephson cur- 
rent from each of the components of the anomalous Green's 
functions. One may rewrite Eq. ( T23T l as 

j = jt + js + Jesp, (25) 

where {jt, j s , j'esp} represent the contribution to the Joseph- 
son current from the S z = triplet, the singlet, and the 
equal spin-pairing anomalous Green's functions, respectively. 
These are defined as 

/oo 
d £ Re{2Af t (e)}tanh(/3 £ /2), 
- OO 

/OO 
deRe{2M s (e)} tanh(/?e/2), 
-OO 

/OO 
deRe{M T (e) + M;(e)} tanh(/?e/2), 
-OO 

(26) 

with the definition Mo = NpeD /2 and 

Now, we get d x M t = X {M S - M a ) = {a =f, |) for 
4> = tt/2 from Eq. (17). With Eq. (20), we have M t = and 
similar M s — M a = by virtue of Eq (21) and the fact that 
c l,r{s) = c*l r{—s)- Therefore, the total Josephson current 
j becomes zero for <fi = tt/2. In fact, an equivalent analytical 
approach is viable to show that the current vanishes in the case 
h = 0. In that case, one finds that d x M t = d x M s = and 
M t — at x — and M s = at x = d by means of the 
boundary conditions and the Usadel equation. 

At this point, it is important to underline that al- 
though the magnitude of the Josephson current in an even- 
frequency/F/odd-frequency Josephson junction depends on 
the orientation of the exchange field, our main message is that 
no current can flow in the absence of a field while the pres- 
ence of an exchange field in general induces a current except 




FIG. 7: (Color online) Plot of the current-phase relationship for sev- 
eral values of cj>. We have chosen h/Ao = 30 and d/£ = 2.0. 

for the special case where the exchange field is parallel to the 
spin of the Cooper pair, (f> = tt/2. 

IV. DISCUSSION 

In our calculations, we have neglected the spatial varia- 
tion of the pairing potential near the interfaces. This is per- 
missable for either low-transmission interfaces or if the super- 
conducting region is much less disordered than the F. 2 Also, 
we have considered non-magnetic interfaces, which are rou- 
tinely used in experiments. Including spin-flip scattering in 
the normal region is not expected to alter our qualitative con- 
clusions since spin-flip scattering alone cannot induce triplet 
pairing in a normal metal in proximity to an even-frequency 
superconductor in the diffusive limit 21 . Our choice of study- 
ing the diffusive limit ensures that one may disregard the 
generation of possible odd parity symmetry components of 
the superconducting gap 12 , which could have caused ambi- 
guities in the interpretation of experimental results obtained 
in our proposed setup. We expect that the predicted effect 
could be experimentally observed in disordered superconduc- 
tor/ferromagnet/superconductor junctions using superconduc- 
tors with different symmetries with respect to frequency. The 
junction widths need to be a few coherence lengths, which is 
well within reach with present-day technology. 

V. SUMMARY 

In summary, we have proposed a method of identifying 
highly unusual superconducting states with the conceptually 
and technologically important property that the order param- 
eter is odd under exchange of time-coordinates of the elec- 
trons constituting a Cooper-pair. Remarkably, we find that 
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an exchange field quite generally induces a Josephson effect 
between even- and odd-frequency superconductors. This con- 
stitutes a clear-cut experimental test for such an unusual su- 
perconducting state. Since our qualitative findings rely on 
symmetry consideration alone, they are expected to be quite 
robust. 
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APPENDIX A: JOSEPHSON CURRENT WITH 
SPIN-DEPENDENT SCATTERING 

We here provide some additional details of our calculations 
and also outline how spin-dependent scattering may be taken 
into account in the analytical expressions. We employ the 
linearized Usadel equations under the assumption of a weak 
proximity effect. This assumption is justified in the low- 
transparency regime (tunneling limit), where the depletion of 
the superconducting order parameter near the interface may 
also be disregarded. In the superconducting reservoirs, we 
employ the bulk solution which reads 



ffsc 



/ 
V- 



C 3 




s,e" 



l x., 







-OtjSjB 














(Al) 



where Cj — cosh^-, Sj = sinh^j, 8j = atanh(A :) (e)/£), and 
j = L,R denotes the left and right superconducting region. 
Here, Xj denotes the broken U(l) phase in superconductor j, 
and we use the convention \R — ~Xl = X- F° r an odd- 
frequency superconductor on side j, we have a.j = 1 while 
for an even-frequency superconductor on side j, one has ctj = 
-1. 

The Kupriyanov-Lukichev 15 boundary conditions now read 



2jdgd x g = -[<?, <?l] U=o, 
2-fdgd x g = [g,g R }\ x=d - 



(A2) 



In the normal region, we may write the Green's function as 

I / t R (e)a+/s R (^i2 

9 \-\f?(-e)n + f?(-e>l£* I 

(A3) 

where the subscripts 't' and V denote the triplet and singlet 
part of the anomalous Green's function. Note that the triplet 
part is odd in frequency. In the following, we will consider 



an exchange field h || z, i.e. perpendicular to the spin of the 
Cooper pair, such that there are no equal-spin pairing compo- 
nents f„ (er =|, 1) of the anomalous Green's function. Intro- 
ducing f± = ft ± / s , we may write the boundary conditions 
more explicitely. At x = one obtains 



2jdd x f+ = 2f+c L - 2 SL e" IX , 
2-/dd x f- = 2f_c L - 2a LSL e- IX , 

while the same procedure at x — d yields 

2 1 dd x f+ = 2s R e lx - 2f+c R , 
2jdd x f_ = 2a R s R e 1 * - 2f_c R . 



(A4) 



(A5) 



The linearized Usadel equations in the normal region may be 
formally obtained by assuming that |/±| <C ljS^- 1 



d 2 x f t ±d 2 J s +A±f t ±B±f s = 0, 
where we have introduced 



(A6) 



2i(e ±h) - g so 



g s fS z 



2i(e ±h)- 



gsi{2S xy + S z ) 



(A7) 



It is possible to find a general analytical solution for the func- 
tions {/ t , / s }, and this can be written as 



- — k- j 



ft = cie '— + c 2 e"~ + c 3 c 
_ Cl S-c- k - x +c 2 S-e k - 

fs — 



-k 4 



c 4 e fc + 3 
,-fc 



c 3 S + e- k + x + c A S+e k 



2(£_ - B^ 



(A8) 



where {cj} are constants to be determined from the boundary 
conditions Eq. (lA4t and iA5i . Also, we have defined the 
auxiliary quantities: 



k± = -y/(-A+-A--B+-B-)±R, 

S± = {A + +A_ - B + - B_)± R, 
R = [{A+ + A-) 2 + (B+ + BJf - 4(A+B_ + A_B+) 
+ 2(A + - A-)(B + - B-)] 1 / 2 . (A9) 

First, we note that 

f ± = Cl e~ k - X (l ± L_) + c 2 e k - x (l ± L_) 

+ c 3 e- k + x (l ± L+) + c A e k+x {\ ± L+), (A10) 

with the definition L± = iDS±/(8h). We also introduce 



1 for 
«l for 



(All) 



and similarly for L — > R. After lengthy calculations, we fi- 
nally aiTive at an explicit expression for the coefficients {c^}: 
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G E IE G 

c 4 = J" ~' t , ~ , c 3 = ~(G+ + c 4 F+)/E+, c 2 = Y a + y 3 c 3 + Y 4 c 4 , a = X + X 2 c 2 + X 3 c 3 + X 4 c 4 . (A12) 

We have defined the auxiliary quantities: 

= 1 ~ L - x = 5 L+ s L c- lx = 7 rffc_ - c L 

1 + L-' ° (l + L_)(cL+7dfc_)' 2 7 ^-+cl' 
(l + L + )(c L + 7 dfc + ) _ (l + L + )( 7 rffc + -cl) 

3 (i + i_)( Ci + 7 dfc_)' 4 (i + L_)( 7 dfc_ + CL y 

8 L e- 1 HPd L + - S L -) 



(A13) 



Y 
Y 3 
Y 4 



( Ci -7dfc_)[P(l + L_)-(l-L_)]' 

(1 - £+)(c L + 7 rffc + ) - (1 + Z-+)(Pc L + P 7 (ifc + ) 

(c L -7rffc_)[P(l + L_)-(l- J L_)] 
(1 - L + )(c L - jdk+) - (1 + L+)(Pc L - P 7 ^+) 



(c L -7dfc_)[P(l + L_)-(l-L_)] 

(A14) 

£± = 2c- fe + d (l ± L+)(c K - 7 rffc+) + 2(1 ± i_)[r 3 c fc - d (c K + 7 dfc_) + (X 2 Y 3 + X 3 ) C - k - d (c R - 7 dfc_)], 

P± = 2c fe + d (l ± L+)(c K + 7 rffc+) + 2(1 ± L_)[r 4 c fe - d (c K + 7 dfc_) + (X 2 Y 4 + X 4 ) C - k - d (c R - 7 dfc_)], 

G± = -2^±s fl e IX + 2(1 ± L_)[r c fe - d (c K + 7 dfc_) + (X 2 Y Q + X Q ) C - k - d (c R - 7 dfc_)]. (A15) 

The above equations may be considerably simplified by considering only uniaxial spin-flip scattering. Setting the planar spin-flip 
and spin-orbit scattering rates to zero, we obtain the anomalous Green's function as 



f± = c 1± e— + c 2± e-^, k ± = f l{£± !^-^\ (A16) 

with the coefficients {p a = vydk a , a = ±) 

SLaSLe~ lx e lk " d (p a + c R ) + 5 Ra s R e lx {p,j - cl) c 2<y (p a + c L ) - S Lcr s L e~ lx 

c lk ° d (p a + c L )(p a + c R ) -c lk ' d {p a - c R )(p a - c L ) p a - c L 

Once g R has been obtained, the Josephson current may be calculated according to the formulas in the main text. 



APPENDIX B: PAULI MATRICES 



The Pauli-matrices used in this paper are defined as 

^=(i o)'- = C 1 )'^=(o -l)' 

1= (J !)' i = (oi)'^(o S)' 

h = (a t) ' h = (4 ' h = (o -l) ' 
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